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Abstract 

The evolution of magnetospheric indices on temporal scales shorter than that of substorms is char- 
acterized by bursty, intermittent events that may arise from turbulence intrinsic to the magnetosphere 
or may reflect solar wind-magnetosphere coupling. This leads to a generic problem of distinguishing be- 
tween the features of the system and those of the driver. We quantify scaling properties of short term (up 
to few hours) fluctuations in the geomagnetic indices AL and AU during solar minimum and maximum 
along with the parameter e that is a measure of the solar wind driver. We flnd that self-similar statistics 
provide a good approximation for the observed scaling properties of fluctuations in the geomagnetic in- 
dices, regardless of the solar activity level, and in the e parameter at solar maximum. This self-similarity 
persists for fluctuations on time scales at least up to about 1 — 2 hours. The scaling exponent of AU 
index fluctuations show dependence on the solar cycle and the trend follows that found in the scaling 
of fluctuations in e. The values of their corresponding scaling exponents, however, are always distinct. 
Fluctuations in the AL index are insensitive to the solar cycle as well as being distinct from those in 
the e parameter. This approximate self-similar scaling leads to a Fokker-Planck model which, we show, 
captures the probability density function of fluctuations and provides a stochastic dynamical equation 
(Langevin equation) for time series of the geomagnetic indices. 
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1 Introduction 

O , 

The Earth's magnetosphere can be considered as non-linear, dissipative system which is driven by the 
J>^time varying solar wind. Accumulated energy is ultimately dissipated, at least in part, through a system 
of currents generated in the auroral zones of the ionosphere. These currents produce small changes in 
• • the terrestrial magnetic field which can be used to monitor magnetospheric activity. The complex be- 
havior of the magnetosphere, as suggested by many observations (see, for example, [Norton et aZ.(1999), 
Lewis{199l), Sitnov et aZ.(2000), Takalo et a/.(2000), Vassiliadis et a/.(2000), Voros et al.{2002)]), could 
then be attributed either to intrinsic magnetospheric processes, the complex nature of its coupling with 
' the solar wind and the ionosphere or both. 

Recent observations suggest that the multi-scale nature of this coupling is a fundamental aspect of mag- 
netospheric dynamics (see, for example [Chang{1992) , Chapman and Watkins{2001) , Klimas et a/.(1996), 
Ukhorskiy et a/.(2003), Vassiliadis et a/.(2003), Weigel et a/.(2003)]). Evidence is provided by a variety of 
observations which exhibit statistical properties previously identified as hallmarks of multi-scale systems. 
For example, bursty transport events have been reported in the magnetotail[Angelopoulos et aZ.(1992)] and 
their auroral signatures suggest self-similar statistics [Lui et a/.(2000), Uritsky et a/.(2001)]. The fiuctua- 
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Figure 1: Exponents of unconditioned generalized structure functions as a function of order m for fluctu- 
ations in the e(o), AU{o) and AL{A) index during solar maximum. The inset shows structure functions 
Sm of orders m = 1 — 4 for the e parameter. 



tions in the ground based measurements of the magnetic field are non-Gaussian and also exhibit scaling 
[Consolini et aZ.(1996), Kovdcs et a/.(2001), Voros et aZ.(1998)]. In the context of time series analysis, ge- 
omagnetic indices arc of particular interest as they provide a global measure of magnctospheric output and 
are evenly sampled over a long time interval. These indices also show non-Gaussian statistics of fluctua- 
tions and anomalous scaling over the short time scales of up to few hours [Consolini and De Michelis{1998) , 
Takalo et a/.(1993), Hnat et aZ.(2003a), Stepanova et a/.(2003), Tsurutani et aZ.(1990)]. The extent to which 
observed statistical features of the geomagnetic indices arise directly from those of the solar wind driver or 
the auroral currents is of fundamental interest. This is an example of the generic problem of distinguishing 
between features intrinsic to a driven system and those in the driver, when both show scaling. Some recent 
studies has focused on direct comparison of scaling properties of the driver with these found in the ge- 
omagnetic indices [Freeman et aZ.(2000), Uritsky et aZ.(2001), Hnat et aZ.(2003a)] to establish whether, to 
the lowest order, they are directly related. 

The difficulty with interpreting these observations arises from the fact that statistical features described 
above can be recovered from many existing models. Self-Organized Criticality (SOC) and turbulence have 
both been extensively used [Angelopoulos et a/.(1999), Consolini and De Mic/ie/is(1998)], 
[Kozelov and Kozelova{2{){)'i) ^ Uritsky and Pudovkin[\'d9S)] in the past. Practically, one needs to obtain 
experimental constrains with which different models with similar characteristics can be tested. In this 
paper we present one possible approach to characterizing the time series in the context of scaling that 
docs not rely on a specific model of multi-scale systems [5'ornette(2000), Hnat et a/.(2003a)]. The aim is to 
obtain statistical scaling properties directly from the data. 

Here, we will examine the statistical properties of Akasofu's e[Perreault and Akasofu{191S)] parameter, 
which represents the energy input from the solar wind into the magnetosphere, and that of magnctospheric 
response as seen by the geomagnetic indices. Previously, scaling has been quantified over a 10 year data set 
for the indices and a comparison between 5e and the indices included, but was not restricted to, the solar 
minimum (1984 — 1987) [i/nai et a/.(2003a)]. Here, we will perform this comparison over intervals of solar 
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Figure 2: Structure functions of orders m = 1 — 6 for fluctuations in the e parameter at solar maximum. 

minimum and maximum separately. The statistical description of the experimental data will be extended 
to 10 standard deviations of the fluctuations. 

To facilitate the comparison of all considered quantities we will first explore to what extend their 
fluctuations exhibit approximate self-similar scaling for temporal scales of 1 — 2 hours. The quality of this 
self-similar approximation combined with values of the scaling exponents obtained at the solar minimum 
and maximum can be used to characterize each quantity. We will see that values of scaling exponents 
on these temporal scales for the geomagnetic indices are different from these found in the solar wind e 
regardless of the phase of the solar cycle. Remarkably, the scaling exponent of the AL index is unchanged 
between solar minimum and maximum whereas the AU scaling exponent changes with the solar cycle. In 
this respect, the AU index seems to follow the trend found in the driver, e i.e., the value of scaling exponent 
increases with increasing solar activity. We then construct a Fokker-Planck model for fluctuations in the 
geomagnetic indices and the e at solar maximum as these exhibit the most satisfactory self-similar scaling. 
This allows us to obtain analytically a functional form of the fluctuation probability density function 
(PDF) which we can then check against the data. A stochastic dynamical model can then be formulated by 
considering the most general form of the Langevin equation and deriving functional forms of the coefficients 
that are consistent with the Fokker-Planck equation (see, for example [Hnat et aZ.(2003b)]). 

2 Data and Methods 

2.1 Data sets 

To facilitate this analysis we used multiple data sets that spanned over different phases of the solar cycle. 
Two year intervals of data were selected centered on solar minimum and solar maximum. The solar wind 
data were obtained from WIND and ACE spacecraft observations. These were collected in the vicinity 
of the LI point approximately 1 AU (Astronomical Units) from the Sun. The periods of coverage, final 
sampling frequencies and number of samples are given in the Table 1. The geomagnetic indices and the 
corresponding spacecraft data sets are not contiguous. The calibrated geomagnetic data set, from which 
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Figure 3: Same as figure 2 for the AU index. 

intervals of interest has been selected, spans from January 1978 to December 1988 inclusive, while the 
spacecraft data are available starting from 1995 for WIND and 1998 for ACE. This available data coverage 
does not permit examination of successive solar cycles. We thus need to assume that the statistical 
properties of fluctuations are invariant from one solar cycle to the next. In the case of the spacecraft data, 
these include slow and fast solar wind streams. 

The solar wind velocity measurements, provided by the SWE instrument on board of WIND and ACE 
spacecraft, have varying temporal resolutions. In the case of WIND this resolution is in the range of 
75 — 98 seconds while for the ACE spacecraft it changes between 60 — 120 seconds. The magnetometer data 
sets, on the other hand, have fixed temporal resolution of 46 seconds for WIND MFI instrument and 16 
seconds in the case of the ACE magnetometer. The SWE data sets have been then re-sampled using linear 
interpolation to give uniform resolution of 92 seconds for WIND (twice the magnetometer resolution) and 
64 seconds for the ACE spacecraft (four times magnetometer resolution). No other post processing, such 
as detrending or smoothing, was applied to data. The e parameter is defined [Perreatt/i and Akasofu{1978)] 
in SI units (Watts) as e = v{B^/iJ,o)lQsiii^{@/2), where Iq ^ 7Re and © = a,TCtaiL{\By\/Bz), and was 
calculated from WIND and ACE spacecraft key parameter databases. 

All techniques discussed here are based on differencing of the original time series over a range of 
temporal scales r. This method is often used in turbulence studies to compare the properties of fluctuations 
on different spatio-temporal scales (see, for example, [Frisch{1995)]) . For a given time series x{t) a set 
of differences 6x{t,T) = x{t + r) — x{t) will then capture fluctuations on temporal scale r. Here, we will 
examine the statistical properties of the PDF of fluctuations Sx{t,T). The r parameter will be given in 
power law form such as r = (5t^t/(1.2)"' seconds, where StAU is a sampling time of the AU time scries (here, 
1 minute) and n > 1 is an integer. This choice of r gives a uniform distance between points when plotted 
on the logarithmic scale while the small base of the power law (1.2) assures that the adequate number of 
temporal scales are explored. We stress that the differencing is performed only if both x{t + t) and x{t) 
exist and are separated by time interval r. 
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Figure 4: Same as figure 2 for the AL index. 



2.2 Statistical Methods 

Generalized structure functions (GSF) Sm are widely used to characterize non-Gaussian processes 
[Frisch{l995) , Hnat et oZ.(2003a)]. These functions can be defined for fluctuations 6x{t,T) as 5'm(r) = 
dfel™"), where m can be any real number, not necessarily positive. If Sm exhibits scaling with respect 
to the time lag r we also have r^("^). In this case a log-log plot of Sm versus r should reveal 

a straight line for each m and the gradients correspond to values of Ct'^)- Generally, CC"^) can be a 
non- linear function of order m, however if CI^t-) = Q;m {a constant) then the time series is self-similar 
(or more precisely, self-affine) with single scaling exponent a. This special case leads immediately to a 
Fokker-Planck description [Hnat et aZ. (2003b)]. The difficulty with computing GSF for higher orders, say, 
m > 4 arises from the slow convergence of this method and its sensitivity to large statistical errors in 
extremal events in the tails of the distribution. These effects can, as we shall see, lead to large errors 
in the estimation of C("^) (see also [Horbury and Balogh(1997)] for the discussion of error estimation for 
structure functions). One possible approach is to eliminate these extreme events from the fluctuation 
time series 5x(t, r) in a way that is consistent with the growth of the self-similar fluctuations' range on 
each temporal scale. This process is referred to as conditioning. Previously, a similar technique based on 
the wavelet filters has been used to separate the intermittent parts of the signal from the homogeneous 
noise in the AE index ddXa[Kovdcs et aZ.(2001)]. We will condition our GSFs by imposing a threshold A 
on the fiuctuation sizc[Hnat et aZ.(2003a)]. The threshold will be based on the standard deviation of the 
differenced time series for a given r, A(^t) = 10(7(t). Under conditioning, the GSF can be expressed in 
term of the fluctuation PDF as: 

<\Sxr >= \SxrP{Sx,T)d{Sx). (1) 
J-A 

This procedure is then consistent with scaling ({m) = ma if it is present in the data, but for threshold A 
sufficiently large it does not enforce it on the data. 

The PDF rcscaling technique is a generic and model independent method of testing for statistical self- 
similarity in the data set. If the data is self-similar, then a single argument representation of the fluctuation 
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Figure 5: Exponents of conditioned generalized structure functions as a function of order m for fluctuations 
in the e(o), AU{o) and AL(A) index during the solar maximum. 
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Figure 6: Structure functions Sm of orders m = 1 — 6 for fluctuations in the e parameter at solar minimum. 
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Figure 7: Same as figure 6 for the AU index. 

PDF, Ps{dxs), can be found in the form: 

P(fe,r) =T-"P,(5xr-°), (2) 

where a is the rcscahng exponent. Substituting the rescaled quantities Ps and 6xs = Sxt~°' into the GSF 
definition given by (1) we obtain: 

< \§xr >= r'"" r° \5xsrPiSx,)di6xs) oc r^^"*), (3) 

where the integral now has no exphcit dependence on temporal scale r. This then immediately relates the 
PDF rescaling to self-similar scaling of the GSF with C("^) = "t-ci- 

In this approach PDFs are generated using non-overlapping intervals of the original data, i.e., 6x{t, r) = 
a;[mT] — x[{m — l)r]. The method assures that fiuctuations are not temporally correlated-an important 
assumption for a Fokker-Planck model wc will consider later. These two methods are thus complementary as 
one provides a scaling exponent while the other gives an underlying probability distribution of fiuctuations. 

3 Results and discussion 
3.1 GSF analysis 

We will first present scaling properties of the GSFs for the indices and the e parameter during solar minimum 

and solar maximum. To illustrate the effect of conditioning, we first show, in the inset of figure 1, a log-log 
plot of structure functions Sm obtained for fluctuations in the raw time series of the AU index at solar 
maximum for orders 1 < m < 6. We see that, for orders m > 3 there is no clear evidence of scaling-the 
points do not lie on straight lines. Similar lack of scaling was also found for fluctuations in the AL index 
and the e parameter. The main panel of flgure 1 shows exponents C("^) obtained by performing linear flts 
to logarithms of moments log[Sm{log{T))]. We see that the curves C(™) are not monotonic functions of 
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Figure 8: Same as figure 6 for the AL index. 

m, excluding tfie possibility of multi-fractal scaling. We now condition this data set as discussed above, to 
check if true scaling properties are not obscured by poor statistics of extreme and very rare events. 

Figures 2-4 show a log-log plot of structure functions for the (5e, 6{AU) and 6{AL) indices at solar 
maximum and for order m from 1 to 6. The main indication of successfully recovered scaling after the 
conditioning process is the quality of the linear fit to log[Sm{T)] versus logir). We clearly recover a family 
of straight lines with slopes C(^)) up to order m = 6 for flTictuations in AU and AL indices and in e. 
This scaling extends up to temporal scales of ~ 1 to ~ 2 hours in good agreement with these reported 
earlier [TaA;aZo et a/.(1993), Tsurutani et aZ.(1990), Hnat et aZ.(2003a)]. Figure 5 shows that fluctuations in 
all quantities, at solar maximum, exhibit approximate self-similar scaling to within statistical errors. The 
size of these error bars combined with the convex shape of the function Cl?^) for the AL index also allows 
a weakly multi-fractal interpretation of the scaling. We stress, however, that the error bars in figure 5 
do not include many other uncertainties (not statistical) that are difficult to estimate. For example, the 
WIND spacecraft magnetometer data has absolute accuracy of about O.lnT and the indices data have 
integer values (also in units of nT). Such discreteness in the time series may lead to erroneous estimates 
of low order moments while the finite size of the data could alter true scaling of the high order moments. 
Independent of any given choice of a model for the functional form of ({m) we can perform a direct 
comparison between the Ci''^) measured for the different quantities at solar minimum and maximum. In 
order to develop a Fokkcr-Planck approach we will then make a further step and assume that a reasonable 
approximation is given by C("^) = "^ct, that is, self-similar scaling. 

Figures 6-8 show structure functions Sm for all quantities at solar minimum and with order m varying 
again from 1 to 6. We see that moments of fluctuations for the geomagnetic indices show satisfactory scaling 
up to temporal scales of 1 — 2 hours. In the case of e at solar minimum there is a departure from a single 
set of scaling exponents C(^) for the smallest time scales r < 12 minutes. To facilitate a comparison with 
conditions at solar maximum and with the indices we will fit straight lines to obtain ({m) for r = [12, 90] 
minutes bearing in mind that this does not capture the behavior of fluctuations on the smallest time scales. 
This change in scaling properties for e may reflect differences between solar wind evolution at solar minimum 
and maximum related to physical properties of slow and fast wind components[Pa5e/ and Balogh{2001)]. 
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Figure 9: Same format as figure 5 for data around solar minimum. 

Figure 9 is constructed identically to figure 5 and shows scaling exponents ({m) for solar minimum. 

To make a comparison between behavior at maximum and minimum we plot, in figure 10, exponents 
C{m) at solar minimum and maximum overplotted for AL, AU and e respectively. Examining these figures 
we conclude that the scaling properties of the AL index fiuctuations are remarkably insensitive to the 
change in solar activity. The values of C(^) ^-^id the corresponding scaling exponents are the same, to 
within the statistical error for solar minimum and maximum. On the other hand the scaling exponents of 
fluctuations in both e and the AU index vary with the solar cycle. The scaling of 5{AU) is distinct from 
these of 6{AL) and Se but follows the trend of 6e. A possible interpretation of these observations is that the 
AL index fiuctuations more closely refiect the internal dynamics of the magnetotail and arc insulated from 
solar cycle related changes in the solar wind. In contrast, the AU index is more strongly coupled to solar 
cycle associated changes in the solar wind driver. This is consistent with our understanding of the global 
roles of these indices (eg., [Baumjohann and Treumann{1996)]) . The fluctuations in AU, however, have 
values of scaling exponents different from that observed for the driver e at solar minimum and maximum, 
which may suggest that (i) e does not completely capture all relevant information about the driver, (ii) 
the indices do not fully capture the magnetospheric response or (iii) the difference reflects the non-linear 
nature of the solar wind-magnetosphere coupling. 

If we compare the scaling exponents of 6{AU) and Se during solar minimum and maximum we see that 
both quantities follow a similar trend. The exponents have values closer to that of Brownian motion (0.5) 
during solar maximum as compared to minimum. Closer examination of scaling exponents for fluctuations 
in the e and the AU index reveals that the difference 0.06 is almost identical for solar minimum 

and solar maximum period, to within the statistical error. This could indicate that the "conversion rate" 
of fluctuations in the driver to those in the AU index is nearly constant and independent of the strength 
of the driver. 

All scaling exponents a derived by fltting C{m) = ma are given in the Table 2 together with the 
approximate maximum temporal scale Tmax for which self-similarity can be identified in the differenced 
time series. These temporal scales have been derived using goodness of fit analysis for moment with 
m = 2. We have also verified that GSF analysis of combined time series over solar minimum and maximum 
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Figure 10: Comparison of functions C{m) during solar minimum and maximum for fluctuations in all three 
quantities: {a)-AL index, {h)-AU index and (c)-e. 

recovers results presented in our previous woTk[Hnat et a/. (2003a)]. 
3.2 Probability Density Function (PDF) Rescaling 

We now present the results of the PDF rescaling analysis which allows us to compare directly the PDFs 
of the studied parameters at solar minimum and maximum. Due to the rather poor scaling of the higher 
moments for the e fluctuations at solar minimum we will not apply this rescaling to their PDFs. We simply 
state that our previous work suggested that PDFs of fluctuations in the geomagnetic indices and that of 
the e fluctuations differed significantly when considered time interval spanned more then a solar minimum 
[Hnat et (2003a)]. 

Figures 11 and 12 show rescaled PDFs of fluctuations in the AU, AL index respectively for solar 
minimum (empty symbols) and solar maximum (fllled symbols) while flgure 13 shows these PDFs for the 

e parameter but only at solar maximum. These PDFs correspond to the function Ps{Sxs) in equation 
(2). Each plot shows overplotted Ps{Sxs) at four temporal scales, r = 10, 16, 26 and 42 minutes. These 
figures show data up to 10 standard deviation on any given temporal scale-consistent with the conditioning 
procedure described above. All rescaling exponents a used to construct these plots, are taken directly from 
the GSF analysis. We find that, when solar minimum and maximum data sets are taken separately, these 
PDFs collapse on a single curve after rescaling (2) is applied. The quality of the collapse for the PDFs was 
checked using the Smirnov-Kolmogorov [Press et a/.(1988)] test and the significance of the null hypothesis 
(both curves drawn from the same distribution) was always found to be above the 0.975 =b 0.05 level. 

The rescaling confirms what we have already found by applying GSF analysis, in that a single exponent 
a is sufficient to give close correspondence of the curves. As we have shown in equation (3) this is consistent 
with approximate self-similar scaling C(^) = "^ck from the GSF analysis. Once rescaled, using the values 
of exponents obtained separately for solar minimum and maximum, we see that the curves are distinct and 
the difference is most clear for the AU index shown in figure 11. 

We also compared the functional form of these curves by applying normalization to their respective 
standard deviation on a given temporal scale, (Ts(t). We found that the normalized PDFs for maximum 
and minimum are indistinguishable within the errors for AU and AL. Similar results were reported for 
ground based measurements of the magnetic field [Weigel and Baker{2003)] where authors also found that 
the statistics of fluctuations, when normalized to the standard deviation, is not sensitive to changing solar 
wind conditions. 
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Figure 11: Rescaled PDFs of the 6{AU) during solar minimum (empty symbols) and maximum (filled 
symbols). Symbols correspond to temporal scales r = 10, 16, 26 and 42 minutes. The dashed line 
represents a solution of the Fokker-Planck model (6) with parameters given in Table 3. 




Figure 12: Same format as figure 11 for the AL index fluctuations. 



11 



4 The Fokker-Planck Model 



The Fokkcr-Planck (F-P) equation provides an important link between statistical properties of the system 
and the dynamical approach expressed by the Langevin equation[iian Kampen{1992)]. The F-P approach 
can be readily applied if fluctuations are self-similar and statistically independent (uncorrelated) 
[van Kampen{1992)]. The above analysis suggests that self-similar scaling is a reasonable approximation 
to the data. The independent nature of increments is enforced by considering non-overlapping intervals 
for differencing, as discussed in Section II B. 

In the most general form the F-P equation can be written as: 

^ = VsAA{^x)P + B{5x)VsxP), (4) 

where P = P{5x, r) is a PDF for the differenced quantity 5x that varies with time r and A{5x) and B{6x) 
are transport coefficients which vary with 5x. It can be shown that, under the assumption of power law 
scaling A{5x) oc Sx^^'^/" and B{5x) oc Jx^^-*^/", a class of self-similar solutions of (4) can be found that 
also satisfies the rescaling relation {2)[Hnat et a/. (2003b)]. These assumptions combined with the use of 
rescaled variables Sxg = Sxt*^ and Pg lead to the following equation: 

-(<^^«):y^ + + -i^^s)^Ps = c, (5) 

ao d{dXs) ao 

where ao,bo , C are constants and a is the rescaling exponent derived, for example, from GSF analysis. The 
general solution of (5) is given by the sum of homogeneous and inhomogeneous solutions [iJnai et a/. (2003b)]: 

PsiSxs) = ^TT-^exp (-^{\8xs\f/- 



bo |(5xs|«o/&o y bo 
rS.. e xp(g(fa;)V") 
(fe^)l-«o/^o 

where ko is a constant and H{5xs) is the homogeneous solution: 



Jo 



d{5x',) + koH{5xs), (6) 



The simple model described above assumes that self-similar scaling persists for all 5x. This assumption is 
expected to hold for a physical system for a large but finite range of 5x. In particular, it will break down 
as 5xs —>■ giving a singularity in the solution Pg as Sxg O.This singularity, however, is integrable so 
that Psd{Sxs) is finite. 

We have found that fluctuations in the geomagnetic indices in solar minimum and maximum and these 
in e at solar maximum exhibit self-similar statistics to a reasonable approximation. We will now show that 
the functional form of the PDF obtained from the F-P model (6) is a good approximation for the observed 
rescaled distribution Pg^Sxg) of fluctuations shown in figures 11-13. On figures 11-13 we have overplotted 
solution (6), shown by thick dashed line, with a taken to be that obtained from the GSF analysis. Table 
3 gives values of all parameters assumed for each of the plotted solutions. We see that all PDFs shown in 
figures 11-13 are well approximated by their F-P solutions. In the case of the geomagnetic indices some 
departures of the predicted curves from the observed distributions do occur and can be attributed to the 
asymmetry of these observed PDFs. 
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Figure 13: Same format as figure 11 for the e parameter fluctuations at solar maximum. 



We note an obvious departure of our predicted curves from the measured PDF for the smallest fluctu- 
ations, in all considered cases. This is due to the functional form of (7) where H{6xs) oo when 5xs — >• 
arising from the assumption that the self-similar scaling extends to arbitrary small fluctuations. To model 
this part of the curve, we would need to include the scaling, or lack of thereof, introduced by the uncertainty 
in the measurements. We would expect such processes to be dominant for the smallest fluctuations. For 
example, if we assume that the smallest fluctuations are dominated by Normally distributed noise, then a 
diffusion model with a constant diffusion coefficient Dq could, in principle, be used to tame this singular 
behavior. This stochastic approach can be extended to obtain the Langevin equation for the dynamics of 
the fluctuations [i^nai et aZ. (2003b)]. The Langevin equation can be written in the most general form as: 



di6x) 
dt 



/3(fe)+7(fe)e(t), 



(8) 



where the random variable ^(t) is assumed to be ^-correlated. Equation 8 can be transform into purely 
additive noise form: 

^ = M + £(t) (9) 
dt 7(z)+^^^' 



where z = J^^ l/^{5x')d{5x'). It has been shown [Hnat et aZ.(2003b)] that one can obtain a functional 
form of coefficients P{6x) and 'y{5x) in terms of oq, bo (from equation 5) and the scaling exponent a. Such 
an equation provides a dynamical model for time series with the required statistical properties. 



5 Summary 

The response of the Earth's magnetosphere to the solar cycle and, by implication, a changing character of 
solar wind activity, illuminates the interplay between intrinsic magnetospheric dynamics and solar wind- 
magnetosphere coupling. Statistical studies provide a simple and yet unifying way to quantify this behavior 
in the context of models for intermittency. In this paper we considered scaling properties of the solar wind 
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driver, quantified by the e parameter, and geomagnetic indices during solar minimum and maximum. We 
find that: 

1. Fluctuations in the geomagnetic indices show approximate statistical self-similarity for a range of 
temporal scales. Fluctuations in the e at solar minimum show departure from scaling for r <~ 10 
minutes. The self-similar scaling emerges as a reasonable approximation for fluctuations 5e at solar 
maximum. Fluctuations in the geomagnetic indices exhibit self-similar scaling on temporal scales 
between ~ 2 minutes to ~ 1 — 2 hours. The fluctuations in e scales from ~ 2 minutes to ^ 1.5 hours, 
but only at solar maximum. 

2. Fluctuations in the AL index exhibit scaling properties insensitive to the phase of the solar cycle. 

3. The scaling exponent of 5{AU) changes with the solar cycle and the trend follows that of the e 
parameter 

4. The value of the scaling exponents of indices and that of the e parameter differ from each other at 
both solar minimum and maximum. This difference between scaling exponents of 6{AU) and the 
driver Se is approximately the same at solar minimum and maximum. 

5. A Fokker-Planck approach can be used to model the fluctuation PDF for the geomagnetic indices in 
both phases of the solar cycle and the e at solar maximum to a good approximation 

The approximate statistical self-similarity found for the indices for solar minimum and maximum and the 
e at solar maximum is consistent with complex multi-scale processes such as turbulence or Self-Organized 
Criticality (SOC). The distinct values found for scaling exponents may reflect physical differences in the 
solar wind and the magnetosphere but may also be due to the very different way in which these quantities 
are derived. The fluctuations in the AU index depend on the solar cycle but the scaling exponent is distinct 
from that of e fluctuations. Interestingly, the difference between scaling exponents of S{AU) and the driver 
5e appears to be approximately constant. These observations, when combined together, suggest that the 
process involved in generating fluctuations in the AU index is coupled to the solar wind driver, as seen in 
the solar cycle dependence. In contrast to the AU index fluctuations, these in the AL index are nearly 
insensitive to the change in solar cycle implying that the AL index is a measure of activity intrinsic to the 
magnetosphere. This is consistent with the AU index more closely monitoring activity on the day-side and 
AL reflecting activity in the magnctotail. 

The self-similar scaling of fluctuations allows us to model their statistics using a Fokker-Planck ap- 
proach. We obtained analytically a functional form of the fluctuation PDF which approximates the mea- 
sured PDF rather well. We stress that such an approach links the statistical features discussed here to 
dynamical modeling of the time series via stochastic Langevin equations. 
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Table 1: Data sets description 



Quantity 


Solar cycle 


dt[s] 


Dates 


N[mln] 


Source 


AL, AU 


Min 


60 


01/85 - 12/86 


1.05 


WDC STP 


AL, AU 


Max 


60 


01/79 - 12/80 


1.05 


WDC STP 


e 


Min 


92 


08/95 - 07/97 


0.63 


WIND 


e 


Max 


64 


01/00 - 12/01 


0.68 


ACE 



Table 2: Scaling indices derived from GSF analysis 



Quantity Solar cycle a from GSF Tmax[hr] 



6AU 


Min 


-0.35 ±0.01 


- 1 


SAU 


Max 


-0.43 ± 0.01 


~ 1 


SAL 


Min 


-0.39 ± 0.02 


~ 2 


SAL 


Max 


-0.37 ±0.03 


~ 2 


Se 


Min 


-0.26 ± 0.02 


- 2 


Se 


Max 


-0.32 ±0.02 


~ 2 



Table 3: Valuc;s oi paraiiic;lc}rH usc^d for F-P model holulioiis ploltcxl in ligures 11- 



Quantity 


Solar cycle 


6o 


bo/ao 


ko 


C [xlO-^] 


SAU 


Min 


170 


1.875 


0.28 


32.5 


SAU 


Max 


16 


1.875 


0.20 


26.4 


SAL 


Min 


2200 


1.8 


0.36 


7.77 


SAL 


Max 


1000 


1.8 


0.32 


6.66 


Se 


Max 


4 X 10^2 


2.15 


5.3 X 10^0 


2.84 X 10-1° 
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